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Abstract. The numerical renormalization group is employed to study a double 
quantum (DQD) dot system consisting of two equivalent single-level dots, each 
coupled to its own lead and with a mutual capacitive coupling embodied in 
an interdot interaction U', in addition to the intradot Coulomb interaction U. 
We focus on the regime with two electrons on the DQD, and the evolution 
of the system on increasing U' /U. The spin-Kondo effect arising for U' = 
{SU{2) X 5(7(2)) is found to persist robustly with increasing U'/U, before a 
rapid but continuous crossover to (a) the SU{A) point U' = U where charge and 
spin degrees of freedom are entangled and the Kondo scale strongly enhanced; 
and then (b) a charge-Kondo state, in which a charge-pseudospin is quenched 
on coupling to the leads/conduction channels. A quantum phase transition of 
Kosterlitz-Thouless type then occurs from this Fermi liquid, strong coupling (SC) 
phase, to a broken symmetry, non-Fermi liquid charge ordered (CO) phase at 
a critical U^. Our emphasis in this paper is on the structure, stability and 
flows between the underlying RG fixed points, on the overall phase diagram 
in the {U, C/')-plane and evolution of the characteristic low-energy Kondo scale 
inherent to the SC phase; and on static physical properties such as spin- and 
charge-susceptibilities (staggered and uniform), including universality and scaling 
behaviour in the strongly correlated regime. Some exact results for associated 
Wilson ratios are also obtained. 
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1. Introduction 

The advent of semiconducting quantum dots has stimulated in recent years a 
remarkable renewal of interest in the Kondo effect [1], a central paradigm of condensed 
matter science more traditionally associated with bulk systems [2]. At its simplest 
the classic spin- Kondo effect arises [3, 4, 5], a single spin in an odd-electron dot 
being quenched by coupling to the low-energy degrees of freedom of the conduction 
electrons in a metallic lead; as manifest strikingly in e.g. the unitarity limit for the zero- 
bias conductance at low temperatures [5]. There is naturally considerable interest in 
coupled quantum dots at present, notably double dot (DQD) systems which are under 
active investigation both theoretically [6-18] and experimentally [19-27]. In addition 
to their possible utility in mesoscopic devices, circuitry and information processing, 
such systems offer the prospect of novel strongly correlated electron states, reflecting 
the inherent relevance of both spin and orbital degrees of freedom. 

Probably the simplest DQD system is that in which two equivalent 
semiconducting dots, each in effect consisting of a single level and separately connected 
to their own leads/conduction channels, are mutually coupled by a capacitive 
interaction embodied in an interdot interaction U' in addition to the usual intradot 
interaction U. The DQD itself can contain up to n = 4 electrons, controlled by varying 
the dot levels (e.g. by application of suitable gate voltages). The problem has been 
studied quite extensively in the n = 1-electron sector [6, 7, 8], where the degenerate 
DQD states {nL,nii) = (1, 0) and (0, 1) dominate and lead to a beautiful example of 
the SU{A) Kondo effect involving spin and orbital degrees of freedom [7, 8]. 

In the present paper by contrast we report an in-depth study of the n = 2-electron 
regime, and its evolution as the ratio U'/U of interdot to intradot interaction strengths 
is increased. We uncover a rich range of physical behaviour, qualitatively quite distinct 
from the n = 1 sector. For U' = the dots are obviously uncoupled, the DQD states 
{riL, ur) = (1,1) dominate and the normal spin-Kondo effect {SU{2) x SU (2)) prevails. 
That essential behaviour is found to persist on increasing U' /U , until close to U'/U = 1 
where wc find a rapid but continuous crossover to the SU{4:) Kondo state occurring 
at U'/U = 1. Hero the six DQD states (1,1), (2,0) and (0,2) are degenerate, the 
spin and charge degrees of freedom are entangled, and the characteristic low-energy 
Kondo scale is markedly enhanced. On increasing U'/U just beyond the 5?7(4) point 
the system then enters a charge-Kondo state, likewise a Fermi liquid state, where a 
charge pscnidospin arising from the DQD states (2, 0) and (0, 2) is Kondo quenched 
on coupling to the conduction channels. 

That Kondo riuenching is however fragile: for only a small increase in U'/U 
the charge-pseudospin tunneling is suppressed, the associated Kondo scale collapses 
rapidly but continuously, and the system undergoes a Kosterlitz-Thouless quantum 
phase transition at a critical U^ from the Fermi liquid phase to a charge ordered 
phase, a broken-symmetry non- Fermi liquid state with In 2 residual entropy reflecting 
the degenerate, unquenched (2,0) and (0,2) DQD conflgurations. Here we study the 
problem using Wilson's numerical renormalization group (NRG) approach [28, 29, 30] 
as the natural method of choice. A preliminary account of some of the results has 
been given in a recent letter [18]. 

The paper is organised as follows. The model is introduced in section 2, and the 
essential steps of the NRG procedure for it are outlined in section 3. RG fixed points, 
their stability and flows between them in different parameter regimes, are discussed in 
section 4, with particular attention given to the charge ordered (CO) fixed point that 
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Figure 1. Schematic representation of tlie double dot system 



controls the CO phase. In section 5 a number of analytical results are obtained by 
considering how the leading corrections to the two stable fixed points, strong coupling 
(SC) and CO, affect flows in their vicinity. For the SC phase in particular we consider 
the static spin susceptibility, uniform and staggered charge susceptibilities, and the 
linear coefficient of specific heat; and obtain a number of exact results for associated 
Wilson ratios. Numerical results are presented in section 6. We begin by establishing 
the phase diagram in the {U, [/')-plane, together with the evolution of the Kondo 
scale in the various regimes of the SC phase. Detailed results for the T-dependence 
of static properties are then given, notably the 'impurity' entropy, spin- and charge- 
susceptibilities. The important issues of scaling and universality in the SC phase are 
likewise considered. Results for (and issues relating to) single-particle dynamics and 
the linear differential conductance, will be given in a subsequent paper [31]. 

2. Model 

We model the capacitively coupled double dot system by a pair of correlated Anderson 
impurities, connected additionally by an interimpurity Coulomb interaction U' as 
shown schematically in figure 1. In standard notation, the Hamiltonian can be written 
as H = Ho + Hini where 

Ho = J2 ^fc^^L^fc^a + J2en,+ ^(aLc^a + h-c), (2.1a) 

i,k,cr i i,k,a 

Hint = ^Ufii^hii + U'nLfiR (2.16) 

i 

and fii = J2a f'-icr = X^o- cl^c^a number operator for dot i G {L, R} (referring to 

the 'left' and 'right' leads/channels). The first two terms of Hq describe the isolated 
lead and dot levels respectively, while the third term contains the one-electron hopping 
processes between each dot and its associated lead. iJjnt contains the interaction 
terms that render the problem non-trivial; here we include both an intradot Coulomb 
repulsion U between the <j =1- and J,-spin electrons on each dot, and the interdot U' 
that embodies the capacitive coupling of the L and R channels. 

When U' = 0, the Hamiltonian (2.1a,6) is simply a doubled version of the single- 
impurity Anderson model, the spin-Kondo physics of which is now well-understood [2] . 
As we show in due course, the additional interdot correlations introduced by U' > 
tend to disfavour this spin-Kondo behaviour, leading to more complex states in which 
both spin and orbital degrees-of- freedom can be important. It is for example easy to 
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show that when U' = U the model maps onto the 5J7(4) Anderson model [2, 7] 

^ en™ + ^ V{al^c^ + h.c.) + hmfim' (2.2) 

fe,m m k,m m,m' 

with spin and channel indices combined into a single quantity m = {i,<7). The 
behaviour at this SU{4:) point is certainly quite different from the SU{2) spin-Kondo 
physics oiU' = Q (albeit that the associated stable fixed points are fundamentally 
of the same type in this case), and one would like to understand how the system 
evolves between the two with increasing U' . The situation beyond U' = U \s also of 
considerable interest: here one naturally expects charge degrees-of- freedom to become 
more important than spin, as indeed will be seen vividly throughout this paper. 

To obtain a rough idea of its parameter regimes, it is first convenient to look at 
the model in the 'atomic' limit {V = 0, where the dots are disconnected from the 
leads). Here the isolated dot Hamiltonian is 

Hd= {ehj + Uhj-^hii) + U'hLhR (2.3a) 

i=L,R 

= eh+\U'h{h-l) + [U -U'] ^ h.^fm {2M) 

z=L,R 

with n = Hi the total DQD mmibcr operator. For fixed interactions U and U' , 
consider progressively lowering the energies of the dot levels e < (via e.g. a suitable 
gate voltage). With increasing |e|, the total ground-state occupancy of the dots 
takes the form of the familiar Coulomb-blockade 'stairc;ase', the occupancy increasing 
stepwise from n = 1 to n = 4. The points of discontinuity in the staircase correspond 
of course to degeneracy of the n- and (n + l)-occupied configurations, and are easily 
shown to lie at |e| = min(J7, U'), U' + max(f7, U') and U + 2U' . The n = 1 sector 
of the full DQD Hamiltonian (2. la, 6) has been studied previously [8]. In this paper 
by contrast we focus on the n = 2 regime, min(L'', U') < \e\ < U' + max([/, U'), and 
the rich range of physical behaviour arising therein. In particular we consider the 
midpoint, |e| = + U', at which the model is particle-hole symmetric; stressing at 
the outset that the p/i-symmetric point is generically representative of the n = 2 regime 
(movement away from the symmetric point simply generates potential scattering in 
the two leads, of both equal sign and magnitude, and does not alter the essential 
physics of the problem). 

With \e\ = + U' henceforth, there remains of course the whole of the (U, U') 
plane to explore. In the atomic limit, this is clearly divided by the line U' = U. For any 
non-zero Z7, the two electrons occupy different dots in the ground-state when U' < U, 
but share the same dot when U' > U; i.e. labelling the states by their individual dot 
occupation numbers rij, the ground state is {nL,nR) = (1, 1) for U' < U, but (2,0) 
or (0,2) for U' > U, as follows directly from (2.36). A quantum phase transition, 
albeit of a rather trivial kind, thus occurs in the atomic limit. But docs a non-trivial 
quantum phase transition arise when the dots are connected to their leads? In the 
remainder of the paper we answer this, and associated questions, by analysing the 
model within the powerful framework of the NRG [28, 29, 30] . 

3. The Numerical Renormalization Group 

The NRG procedure for the DQD system is naturally quite similar to that of the 
single-impurity Anderson model. We thus refer the reader to [29] for further details, 
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but outline the essential steps here in order to point out the differences that arise 
when the impurities are connected to two conduction bands. 



3.1. Linear Chain Hamiltonian 

Following section II of [29] , this transformation of the Hamiltonian can be divided into 
three main steps. 

(i) One defines a new set of conduction-band operators {(ikia}, in terms of which 
the Hamiltonian takes a one- dimensional, continuum form. Each of the leads is 
assumed to consist of a single flat band, with density-of-states p and bandwidth 
2D. The natural lead-dot hybridization parameter is then the quantity F = 'KpV'^, 
and the DQD Hamiltonian becomes 



1 U' 
+ ^^{£ni + Uhi^hii) + —fiLfiR. (3.1) 

i 

(ii) Both conduction bands are then divided symmetrically about k = into 
logarithmic intervals, the nth of which (for fc > 0) spans the range A~'^"+^^ < 
k < A~". Only a single conduction-band state from each interval — the fully 
symmetric linear combination — couples to the impurity. All other conduction 
band states are neglected at this stage. 

(iii) Finally, one performs a unitary transformation of the conduction-band states, to 
write each lead as a simple linear chain. The transformation is chosen such that 
only one site from each chain couples to the corresponding dot, by defining the 
annihilation operator for this site as 

foia = ^j\dka,,^. (3.2) 

The remaining sites in each chain are described by operators /^^^ (with somewhat 
more complicated expansions in terms of the a|.^^s, see [29]), and the Hamiltonian 
thus takes the form 

1/2 



H / 9V\ ' 

- = ^ i(l + A-') E A-"/'Cn(/L/n+M. + h.C.) + — {fLci. + h.C.) 

j,cr L n=0 ^ ^ 

1 U' 

+ ^ ^{^ni + UhifUii) + —fiLfiR (3.3) 

i 

where ^„ = (1 - A-"-i)/[(l - A-2"-i)(l - A-2"-3)]i/2 (^jth ^„ ^ 1 for n » 1). 



3.2. Iterative Diagonalization 

As originally discussed in [28] , an accurate description of the ground state cannot be 
obtained simply by diagonalising a truncated version of the linear-chain Hamiltonian 
(3.3). This is because the most important conduction band states at low-temperatures 
(those near the Fermi level) map onto the f^^^ sites with the largest values of n. To 
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determine the true nature of the ground state, it is necessary to retain all conduction- 
band chain sites; this is achieved in practice by using an iterative procedure to 
diagonalize the Hamiltonian. 

The first step is to write (3.3) as 



H= Urn i(l + A-i)£)A-(^-i)/2^jv, 



JV— >(X) 



(3.4) 



where the N-site Hamiltonian i?jv is itself given by 



F;V = A(^-l)/2^ 



■JV-1 

E 

n=0 



A-"/'€„(/L/n+l,ia + h-C.) + f l/2(/t.^c,, + h.C.) 



+A(Af-l)/2) 



^(enj + Uhi-^hii) + U'riLriR 



(3.5) 



with dimensionless coupling constants 

^ 2r 



r = 
u = 



1 + A-i 

U' 

1 + A-i/D 



1 + A-iy D 
2 \U' 

'd' 



1 + A- 



The key point here is that satisfies the recursion relation 

F;V+1 = AV2/f^ + Yifl^J^,,,^ + h.C.) ^ T{Hm). 



(3.6) 



(3.7) 



Equation (3.7) has the status of an RG transformation, and lies at the heart of the 
approach [28, 29, 30]. This recursion relation could in principle be implemented 
directly, by diagonalizing Hj^j^i in the direct-product basis of the states of iJjv and 
the 16 states obtained by placing electrons in the {N + l)th orbitals of the two chains. 
In practice, however, the computation time required for the diagonalization of H^^i 
is greatly reduced by making use of its symmetries. As the coupling between the dots 
is purely capacitivc, it is readily shown that Hm commutes with the individual charge 
operators for the L/R subsystems 



Q 



Ni 




p. f 



nicT-' ni<T 



1 



1 



and also the left- and right-hand total spin operators 



N 



n— a, a' 



a.<7' 



(3.8) 



(3.9) 



where a denotes the Pauli matrices. It is thus possible to label each cigcustate of ifjv 
using seven quantum numbers, i.e. IQl, S'l, S'lz; Q;?, 5*^, S'fix; r); whore Si and Siz 
represent respectively the magnitude and ^-component of spin on side i {= L/R), and 
the index r labels the different states of each charge and spin subspace. By using the 
states of i?7v to form basis states for H^^i which are also cigcnstates of QN+i,i, SN+i,i 
and Sz,N+i,i, it can be shown that the H^+i can be diagonalized independently within 
each subspace, thereby reducing the computation time significantly. Furthermore, as 
the states within each spin multiplct arc degenerate, one can use the Wigner-Eckart 
theorem [29] to avoid having to consider each of them explicitly. It is then unnecessary 
to retain the Siz labels, which we thus omit from here onward. 
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It is of course impossible to perforin the iterative diagonalization exactly — 
the dimensionaUty of the Hilbcrt space grows by a factor of 16 at each step. Even 
exploiting the symmetries just described, the complete numerical diagonalization of 
Hn becomes unfeasible after just a few iterations. The solution [28] is to retain only a 
certain number Ns of the lowest cniergy states of _ff,Y to calculate the states of H^+i, 
with A^s chosen as a compromise between reducing the computation time while keeping 
the resulting truncation errors to a minimum. The results in this paper have typically 
been obtained using Ns ^ 20,000 (which we add excludes spin multiplicities). In 
addition, one must also take into account the error resulting from the logarithmic 
discretization of the Hamiltonian. The discretized Hamiltonian (3.3) is an exact 
transformation of (3.1) in the limit A ^ 1. While one might naively imagine that 
A should be chosen as close to 1 as possible, this is not so [28, 29] because the number 
of iterations needed to reach the low-energy conduction band states that dominate low- 
energy properties would become very large; not only increasing the calculation time, 
but more importantly allowing the cumulative effect of the truncation error from each 
iteration to become unacceptably large. It is in fact possible to make corrections for 
the discretization error resulting from even a relatively large A [28, 29]; and for the 
present calculations we take A = 3. 

3.3. Thermodynamics 

The temperature dependence of relevant thermodynamic functions can be calculated 
via the NRG as described in [28, 29]. In this paper, we give results for three 'impurity' 
susceptibilities (borrowing the established terminology associated with Anderson- 
type impurity models). These are the spin, uniform charge and staggered charge 
susceptibilities, denoted by x.,, xt Xc respectively, and defined as 

k^Txs = {gi^nf ( [S, - {S,)]^) (3.10) 

\ / imp 

kBTxf = \([Nl±Nr- {Nl ± Nr)^) (3.11) 

\ / imp 

where Sz is the ^-component of the total spin, and Nl (Nr) is the total charge operator 
for the left (right) lead and dot subsystem (as usual, (O)iinp = (fi) — {Ct)o with (r2)o 
denoting a thermal average in the absence of the impurities/dots). Note that the 
staggered charge susceptibility Xc effectively a 'charge pseudospin' analogue of the 
spin susceptibility Xs, and that at the point U' = U where the model has SU{4) 
symmetry (see (2.2)), Xc = Xs (as discussed further in §5,6). In addition, we shall 
use the impurity free energy F = — kBT In (Z/Zq) (with Zq the partition function in 
the absence of the two dots) to calculate the impurity entropy S = —dF/dT, and the 
linear specific heat coefficient 7 given by C(T) = — ~ jT + C(T^). From here 

on, we will take fcs = g/J-B = 1 for convenience. 

4. Fixed points 

Before discussing the full numerical results, we examine the fixed points that control 
the NRG flows. In particular we shall see that there are two distinct classes of stable 
fixed points, the mere existence of which implies a quantum phase transition in the 
model above a critical interdot interaction U^. 

Two initial remarks on the general nature of NRG fixed points should be made. 
First, the energies of the Hamiltonian (3.5) are well known [28] to depend 
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Table 1. Fixed points of the DQD system, obtained from (3.5) 
by setting T, U and U' to values of either or oo. 



Name 


Abbreviation 


r 


u 


U' 


Free Orbital 


FO 











ST/ (2) Local Moment 


LM-SC/(2) 





oo 





ST/ (4) Local Moment 







oo 


oo 


Charge Ordered 


CO 








oo 


Strong Coupling 


sc 


oo 









fundamentally on whether N is even or odd. One works therefore with fixed points 
of T^; the fixed point Hamiltonian H* satisfying T'^{H*) = H* . The second, more 
subtle point is that each successive RG transformation increases the dimensionality of 
the Hilbert space of the problem; no Hamiltonian can thus remain entirely unaffected 
by T^. But there do exist Hamiltonians, H'f^, for which the low-lying states (i.e. 
the Ith lowest excitations, with I <C N) remain asymptotically unaffected by the RG 
transformation; it is these which constitute the fixed points of the NRG procedure. 

The set of NRG fixed points consistent with the symmetries of the DQD 
Hamiltonian (3.3) may be deduced by setting each of its bare parameters to a value of 
either zero or infinity. This leaves in all cases a pair of 'free' conduction chains, plus a 
piece that describes any remaining dot degrees-of-freedom. Since the free conduction 
chain energy levels themselves converge rapidly [28] with increasing N, while the 
isolated dot degrees-of-freedom are entirely unaffected by the NRG transformation, 
this decoupled Hamiltonian is necessarily a fixed point of T^. 

Table 1 lists the five distinct fixed points resulting from the procedure sketched 
above. The three fixed point Hamiltonians with [/' = are obviously just 'doubled' 
versions of the corresponding Hamiltonians for the symmetric Anderson model [29] 
(although the physical behaviour near these fixed points will generally differ from 
that of the single- impurity case because corrections to the fixed points can include 
terms that couple the two channels). The other two fixed points, LM'^'^^^^ and CO, 
have no analogues in the AIM, and thus generate the distinct physics of the DQD 
model. 

The stability of each fixed point in table 1 can be investigated by linearising 
the NRG transformation in its vicinity. Any perturbation 5H of the fixed point 
that is consistent with its underlying symmetries, is known [28] to have the following 
expansion in terms of the eigenvectors O* , and corresponding eigenvalues A;*, of the 
matrix C* that describes the transformation: 

AiN-i)/2sH = ^ CiXf/^Ot, (4.1) 

(with Ci a constant independent of A'^). By analysing the possible SHa in (4.1), one 
can deduce the eigenvalues A;* by equating the A^-dependences of both sides. If any of 
these SHs generate relevant eigenvectors, the fixed point is unstable. 

For each fixed point of the DQD system, table 2 shows the resultant perturbations 
that generate the most relevant eigenvectors of the linearised NRG transformation. 
Most are straightforward to analyse in practice. The free orbital fixed point is readily 
shown to be tmstable, and the strong coupling fixed point to be stable. Each of the 
local moment fixed points is marginal to leading (linear) order, but our numerical 
results in section 6 show they are both unstable under the full NRG transformation. 
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Table 2. Leading corrections to the fixed points of tfie DQD. The quantity Xj 
represents a spin-i on dot i, while a is an anomalous exponent discussed in the 
text. In this table, we adopt the convention that repeated spin indices a are to 
be implicitly summed over unless otherwise stated. 



Fixed point 


Label 


Perturbation 


Eigenvalue 


Stability 


FO 


SHi 




AV2 


Unstable 




SH2 




A 






6H3 


flLflR 


A 






SH4 




1 


Marginally- unstable 




SH5 




1 


Marginally- unstable 


SC 


SHe 


E.(/L/2..+h.c.) 


A-i 


Stable 




SHr 


'^iiflia flier ^ 1) 


A-i 








Yiiifliaflia ~ 1) 


A-i 




CO 


SHg 


— l)/oi(7/oso- 


1 


Stable if a > 








A-2« 





The charge ordered (CO) fixed point by contrast is more subtle; we thus devote the 
next section to consideration of its behaviour. 

4-1- The charge ordered fixed point 

The CO fixed point Hamiltonian can be obtained by setting {/' = oo in (3.5). This 

decouples the dots from the leads, leaving a pair of free conduction bands plus the two 
degenerate dot charge configurations {iiljUr) = (0,2) and (2,0). To avoid excessive 
use of projection operators in the following, it is implicit that the Hilbert space in 
which the CO fixed point Hamiltonian acts contains these dot configurations only. 

As seen in table 2, two corrections to the fixed point govern its stability. We 
begin by discussing the marginal term 

5^9 = E - ^)fLf0ia =H^L- Ur) J2 {fLjoLa ' fLjoR.) , (4-2) 

which describes potential scattering of the conduction electrons. The key point here 
is that although the magnitude of this potential scattering is equal for both the L 
and R conduction bands, the signs arc opposite and correlated with the two charge 
configurations of the dots. It is this correlation that is responsible for much of the 
interesting behaviour of the DQD model. 

As shown originally in [30] , the eff'ect of potential scattering on the conduction 
band chains can be conveniently incorporated into the definition of the fixed point 
itself. Following the notation of [30] we introduce a coupling constant > to specify 
its magnitude such that the appropriate picture is of a line of CO fixed points with 
different values of K — and thus obtain the following CO fixed point Hamiltonian: 

In the Appendix, we point out that a continuum version of this Hamiltonian is 
precisely the 'full' low-energy model for the DQD system when U' U; and by 



riv-i 



n=0 



4.3) 
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comparison of the discrete and continuum forms obtain a relationship between K and 

the bare parameters of the model, which we show to be in excellent agreement with K 
obtained from the full NRG calculations at large U'. As far as the present discussion 
is concerned, the most important result to be drawn from the continuum version of 
(4.3) is that the potential scattering generates equal and opposite phase shifts in the 
two leads, the magnitude of which can be related to K via 



tan 



7r(l - A-^)j^ 
2 In A 



(4.4) 



which result will shortly prove useful. 

Given that 5Hg above is exactly marginal, the stability of each CO fixed point is 
found to be controlled by the term SHio in table 2, namely SHiq = A + with 



n4 



RrrfoRafoLrr^Lrr- i'^-^) 

This operator describes four-electron cotunnelling processes that interconvert the 
(0,2) and (2,0) dot configurations, and thereby switch the signs of the potential 
scattering on each lead. The relevance of 6Hiq is not obvious from simple power 
counting [12, 32]; the action of (4.5) on the eigenstates of qq{K) is very similar 
to the non-perturbative X-ray edge problem [33, 34] and hence rcqiiires a somewhat 
more sophisticated treatment such as bosonization [35] . By considering the long-time 
behaviour of the correlation function {A'^{t)A)o, one can show that the eigenvalue 
associated with SHw goes as A~^", where a is an anomalous exponent that is related 
to the conduction band phase shift 5 via 

„ = -i.(f-i)\ ,«) 

This anomalous exponent arises also for a particular limit of the Ising-coupled Kondo 
impurity model [12], a point which to which we return in section 5.1. 

The stability of each CO fixed point depends on the corresponding sign of a, 
which can be calculated directly from (4.4) and (4.6). Since K S (0,oo) it is easy to 
see that the line of CO fixed points divides at a critical into two classes: stable for 
K < Kc (i.e. a > 0), and unstable for K > (a < 0). The critical point, at which 
a vanishes, thus corresponds to a phase shift given by 

and hence a critical potential scattering strength 



- 2 In A 

Kc = —r. 7— p- tan 

7r(l - A-i) 



2V ^/2 



(4.8) 



(as we also confirm numerically in the Appendix). Note that in the particle- hole 
symmetric limit considered explicitly, be is entirely independent of the bare material 
parameters. That is not however true in general — deviation from particle-hole 
symmetry generates additional potential scattering of the same sign on the two leads. 
This breaks the simple 'equal and opposite' relationship between the L and R phase 
shifts, and leads to a more complicated form for the anomalous exponent a (but 
as mentioned above does not alter the essential physics, which is wholly robust to 
departure from p-h symmetry). 
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4-2. Quantum phase transition and schematic flow diagram 

Our analysis of the NRG fixed points shows two qualitatively distinct classes of stable 
fixed points. SC and CO. This of course implies the existence of a quantum phase 
transition in the system. Points in phase space that flow to the SC fixed point under 
renormalization comprise the 'strong coupling' phase, where all the degrees of freedom 
associated with the dots arc frozen-out at low temperatures, while the set of points 
which flow to the CO fixed point define a 'charge ordered' phase which has a doubly- 
degenerate ground state. Although the phase transition can of course be approached 
in any direction, we will tend to envisage crossing from SC to CO by increasing U' at 
a fixed U and F; the transition then takes place at the critical point U' = U'^. 

It is helpful at this stage to draw a schematic flow diagram showing how the 
various fixed points of the model arc approached under different choices of the bare 
parameters U, U' and F; this will be useful in the next section when we analyse how 
the various unstable fixed points affect the SC phase, and how the CO fixed point in 
particular controls the quantum phase transition. [The structure of the flow diagram 
is of course deduced from the full NRG calculations, as detailed in section 6.] In 
Figure 2 the flow diagram is drawn in the three-dimensional (Tes/D, Ues/D, U[,q/D) 
effective parameter space (we omit the potential scattering generated in the region of 
the CO fixed points, and thus represent schematically the line of CO fixed points by 
a single point). In this work we are primarily concerned with the strongly-correlated 
regime arising for U /T 3> 1, and hence take the bandwidth D to be the largest energy 
scale; all fiows thus begin near the FO fixed point, as shown in the figure. 

When the bare U' = 0, the flow is flrst to the LM'^^(^) fixed point, and then 
to the stable SC [29]. Switching on a finite U' brings the additional fixed points of 
the DQD system into play, as shown by the successively darker fiow lines in Figure 2. 
When U' = U (i.e. the 'middle' flow line in figure 2), the NRG flow remains in the 
Ucii = U'cs plane throughout, heading to the LM'^'^^^^ fixed point before ultimately the 
stable SC. With any deviation from this SU{4) symmetry, the fiows are drawn away 
from the the U^s = U^g plane, crossing over instead to LM'^^^^^ (for U' < U) or CO 
(for U' > U). The physics around the CO fixed point is in fact crucial to the problem: 
we show in section 6.1 that there is a range of U' above U where the CO fixed point 
is unstable, the flows ultimately turning to SC at large N as depicted by the fourth 
flow in figure 2. At a critical however, the CO fixed point becomes stable and the 
quantum phase transition from the SC phase to the CO phase ensues. In section 6 we 
analyse the numerical results to confirm the above picture, and identify the relevant 
energy scales at which the crossovers between the various fixed points occur. 

5. Fixed point analytics 

First, we focus on results that can be obtained analytically by considering how the 
leading corrections to the SC and CO fixed points affect the flows in their vicinities. 
In particular we show that the effective low-energy Hamiltonian of the system near 
criticality can be obtained from flows close to the CO fixed point, and analyse the 
behaviour near SC to obtain the leading T — > behaviour of the SC phase. 
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Figure 2. Schematic flow diagram for the DQD system. We consider the case 
when the bare ?7 2> F, but with U, U' and T all much smaller than the bandwidth 
D. The diagram shows the effect of increasing U' at fixed U and F (light to dark 
flow lines); at a critical the flow tends ultimately to CO instead of SC. (Note 
that in addition to the three effective parameters Ugfi, U^g and Fefi shown, there 
is an effective potential scattering K generated in the vicinity of the CO fixed 
point.) 



5.1. Effective Hamiltonian near the phase boundary 

We begin with the CO fixed point. For U' U'^, there exists a range of N over which 
the NRG Hamiltonians Hn flow very close to the Hnc of CO fixed points (sec figure 2) 
and can therefore be approximated by H]^ qq plus its leading corrections. To make 
this statement quantitative, suppose that the flow is close to the CO fixed points when 
N = Ni and use (4.3) and (4.5) to write the corresponding Hamiltonian as 



H 



Ni,CO 



A(JVi-1)/2 



n=0 



(5.1) 



(with the constants K and ujiq naturally determined by fitting the energy levels of 
(5.1) to those obtained from the full NRG calculations). 

Suppose that we are only interested in the low-temperature properties of the 
system (where by 'low-temperature' wc mean T < T/v^ ^ DA^^'-^^). These are 
calculated from the sequence of Hamiltonians Hj^,^, H^j^^i, . . ., obtained by applying 
successive RG transformations (3.7) to H^^. While we have here obtained H^^ 
from our analysis of the 'full' DQD Hamiltonian (2.1a,fo), it is clear that any bare 
Hamiltonian that leads to the same Hn-^ under RG flow will show precisely the same 
thermodynamics for T < T^i', such a Hamiltonian thus constitutes an eff'ective low- 
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temperature model of the system. The simplest bare Hamiltonian that generates (5.1) 
is of the form 

+ JzTz (4La«fe'L<T - «Laafc'fla)' (5-2) 

where T represents a spin-^ pseudospin defined by the operators = Yla ^^La'^Ra' 
T~ = Yia '^lia'^La ^^'-^ ~ {fiL — nii)/4,\ and the coupUng constants are Jj_ = 
5(1 + A-i)a;io-D and = (1 + k-'^)kD. The reader can verify that (5.2) gives 
the same i7jvi as (5.1) by substituting it for the full DQD Hamiltonian in (2.1a,6) 
and repeating the subsequent steps of section 3. A new is obtained in place 
of (3.5), which is identical to (5.1) when A'' = Ni. Equation (5.1) thus captures the 
essential low-temperature degrees-of-freedom of the system close to the quantum phase 
transition. 

It is however instructive to recast (5.2) in a somewhat different form. As it stands, 
the pseudospin fiips described by occur in combination with electrons hopping 
between the left and right conduction bands. It is expedient to interchange the spin 
and L/R labels (t<-^ L, J,<-> R) so that in this relabelled form the pseudospin flips are 
accompanied by spin flips of the conduction electrons. With this, (5.2) reduces to 

ijfejCr \ i k,k' I 

This effective low-energy Hamiltonian is generically valid as f/' ^ [/^. A Hamiltonian 
of this form has been investigated previously, in the problem of a pair of Ising- 
coupled Kondo impurities [12]. The latter model displays a quantum phase transition 
of Kosterlitz-Thouless [36] (KT) type, and in section 6.1 we show that the phase 
transition in our DQD system is also of KT type. 

5.2. Thermodynamics of the strong coupling phase 

We turn now to the leading T ^ behaviour of the SC phase, which can be obtained 
from the SC fixed point Hamiltonian H-^ ^q and its most relevant corrections given 
in table 2. We write the effective NRG Hamiltonian as 

Hn,sc « H*^,sc + c^6A-(^-^)/25if6 + a;7A-(^-i)/25if7 + 0JsA-^''-'^/'6Hs, (5.4) 

(with the coupling constants We, and determined numerically). It is 

straightforward but lengthy to calculate the T ^ thermodynamics of this effective 
Hamiltonian; one simply uses the results of section 3.3, treating the fixed point 
corrections to leading order in perturbation theory (see e.g. section V of [29]). The final 
results for the impurity susceptibilities and linear specific heat coefficient as A — > 1, 
are found to be 



_i_ 8 z aoai 
^ Dl + A-^ In A 



-ue - {u)7 ± Us ) 



ai In A 



(5.5) 
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Xs = (giJ-B 



2 2 2 apai 
Ul + A-i In A 

,28 2 aQtti 



-LOq + L07 



ai In A 



(5.6) 



where ao and ai are dependent on the discretization parameter A, viz 

ao = [i(l - A-')] ^ ai = [i(l - A-3)] 5 . (5.8) 

Equations (5.5) (5.7) provide a convenient means of calculating accurately the 
low-temperature thermodynamics from the NRG energy levels obtained numerically. 
More importantly however, their dependences on the three parameters uiq, and 
enable us to deduce some exact results. We shall express these in the language of 
Wilson ratios, defined by 

4 {^kef Xs p± Mjikefxt Q^ 

where we add for later use that Tk = 1/7 provides a suitable definition of the Kondo 
scale Tk characteristic of the SC phase. Here Rs is the double dot analogue of the 
usual 'spin' Wilson ratio, and Rf make up a pair of 'charge' Wilson ratios. It is 
straightforward to show from (5.5)-(5.7) that these Wilson ratios can be expressed in 
terms of the coupling constants of the effective Hamiltonian (5.4) as 

i?. = l-^^. i?.^ = l+ ^"^^"^^ f (5.10) 

ai In A We aim A 

and therefore that the three Wilson ratios are related by 

i?+ + i?^ + 2i?^ = 4. (5.11) 

This A-independent result is satisfied throughout the SC phase, and we believe it to 

be exact. We can moreover go further by focusing on three high-symmetry limits of 
the model, where it is possible to determine the precise values of the Wilson ratios. 
We end this section by considering these limits in turn. 



5.2.1. The SU{2) limit. When U' = 0, the left and right channels of the system are 
disconnected and it is clear that the coupling constant of the effective model (5.4) 
vanishes. We then of course recover from (5.10) and (5.11) the well-known results for 
the AIM [2]: the two charge Wilson ratios coincide (i?+ = R~ = Rc), and Rs+Rc = 2 
for all U. In the strong coupling limit, U ^ 1, the system maps onto a Kondo model 
where each dot is strictly singly-occupied; the charge Wilson ratio Rc then vanishes 
and Re = 2 [28, 37]. 

5.2.2. The SU{4:) limit. At the point U' = U, we can make use of the underlying 
SU{4:) symmetry of the Hamiltonian (see (2.2)); for (5.4) to be likewise SU{4:)- 
invariant, it is necessary and sufficient that ug = 2lvt, such that 

Hn,SC « H*^^sC + ^6A-(^-^)/2 E(/L/2m + 

ra 

(5.12) 

m 

with m = (?', a) as before. Using = 2ti)7 in (5.10) gives Rs = R^ for all U' = U m 
the SC phase. In the strong-coupling limit f/ 1 in particular, the total charge on 
the dots is strictly n = 2, i?+ vanishes again and it follows from (5.11) that [38] 

R^=R- = i?+ = when U' = U ^00. (5.13) 
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5.2.3. The U' — > f7^~ limit. In section 5.1, we showed that the low-temperature 
{T <^ U' — U) properties of the system near the phase transition can be obtained 
from a Hamiltonian of the form (5.1). Here we can use one of the symmetries of 
this Hamiltonian to deduce a relation between and cog when U' is close to U'^ in 
the SC phase. As a result of its reduced dot degrees-of-frecdom compared to the 
full DQD Hamiltonian (3.5), the effective Hamiltonian (5.1) is invariant not only 
under a complete particle-hole transformation, but also under a combined particle- 
hole transformation and interchange of the L and R labels in a single spin channel. 
Such a transformation is performed as follows for the f-spin channel 

with the |-spin operators left unchanged. It is straightforward to show that for Hn,sc 
to be similarly invariant under this transformation, = — 2w7, and thus from (5.10) 

that Rs = R^ as U' U^r . Further, as shown in section 6 the Kondo scale Tk 
vanishes as the transition is approached (and thus 7 diverges), while Xs remains finite; 
hence 

R,=R+ = 0, R- =4: when U' U'' (5.15) 
(holding for any U). 

6. Numerical results: thermodynamics 

We now present numerical results obtained from the NRG calculations. While much 
can been deduced from the general RG framework itself, as above, the numerics provide 
accurate and detailed results for thermodynamics and dynamics over a wide range of 
temperature and frequency scales. We first determine the phase boundary between the 
SC to CO phases in the {U, U') plane, before analysing the Kondo scale that controls 
the low-energy behaviour of the SC phase and vanishes as the SC— >C0 transition is 
approached. Typical temperature dependences of the 'impurity' entropy, spin- and 
staggered charge-susceptibilities with increasing U' are then considered, and shown to 
exhibit universal scaling behaviour in the SC phase. 

6.1. Phase diagram and the Kondo scale 

We determine the phase boundary in the {U, U') plane by a simple bisection method, 
classifying points therein as belonging either to the SC or CO phase by analysing their 
corresponding NRG flows. For a particular choice of ?7 = U /{ttT) and U' = (J/'/ttF), 
the fixed point ultimately reached as ^ 00 can be identified by comparing the 
energy levels of the full NRG calculation with those easily obtained from the single- 
particle energies of the fixed point Hamiltonians (5.4) and (4.3). The energy level 
patterns that arise in the two phases are quite distinct, so the identification of the 
ground state is straightforward. To illustrate this, figure 3 shows the A'-dependence 
of the six lowest energy levels in the {Ql, Sl, Qr., Sn) = (1,0, 0, i) subspace of Hn, 
for two characteristic points 'deep' in the SC and CO phases: U = U/ (ttF) = 7, with 
U' = U'/ivrT) = 6 (SC) and U' = 8 (CO). The initial flows in both phases are clearly 
very similar, as expected because both flows begin near the FO fixed point. With 
increasing N however, the SC and CO energy levels begin to depart from each other, 
due to the influences of the diff'erent fixed points on the two flows. For large N, the SC 
and CO phase levels converge to qualitatively difl[erent limits. The pattern of levels 
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Figure 3. Lowest six energy levels of the {Ql, Sl , Qr , Sr) = (1,0,0, ^) subspace 
of -ffjV) a function of N, at characteristic points deep in the SC (C/ = 7, U' = 6) 
and CO {U = 7,U' = 8) phases. 



SC 



CO 



^ ' ' ' ' ' 

2 4 6 8 10 

U' 



Figure 4. Phase diagram in the {U',U) plane. The solid line is the phase 
boundary between the SC and CO phases, and the dotted line is the S{7(4) line 
U' = U. 



shown for the SC phase in figure 3 is in fact the same for all points in the SC phase, 
while in the CO phase the values of the fixed point energies depend on the potential 
scattering K (see (4.3)), which varies with the choice of bare model parameters. 

Figure 4 shows the resultant phase boundary. As a function of U, the critical 
[/' is finite for all [/ > (and moreover exceeds U, the boundary between (1, 1) 
and (2,0)/(0, 2) configurations in the atomic limit). In strong-coupling in particular 
(U 3> 1), we find — U to be exponentially small, taking the form {U^/U — 1) ~ 
cU-^^^ exp{—n'^U / 16). We show below that extension of the SC phase beyond the 
atomic limit boundary reflects the formation of a charge-Kondo state, where the 
doubly-degenerate charge pseudospin that exists freely in the CO phase is quenched 
by coupling to the leads. Simple physical arguments [18] (reprised in the Appendix) 
show that the stability of this charge-Kondo state depends on the magnitide oiU' — U 
relative to the SU (4) Kondo temperature, and hence we now turn to consideration of 
Tk in the SC phase. 

As stated in section 5.2, we choose to define Tk throughout the SC phase as 
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Figure 5. Evolution of tlic Kondo scale Tk across the SC phase. The main figure 
shows lii(Tx/r) versus U'/U for fixed U = 7 (upper set of points) and f7 = 10 
(lower). For U = 7, the inset shows the Kosterlitz-Thouless behaviour close to 
U'^ ~ 7.046: here, ln(Tx/r) versus {U^ — U')~^^'^ approaches the asymptotic form 
of (6.3) with c ~ 12.80 and a ~ 1.38, shown as a dotted line. 



Tk = 1/7- This is most easily calculated using the analytical result (5.7), with the 
value of uq determined from the numerical energy levels at large N. Figure 5 shows 
the typical strong-coupling behaviour of Tk across the SC phase, as a function of U'/JJ 
for two fixed values oi U (=7 and 10). We see that the SC phase divides naturally 
into two distinct regimes, namely U' < U and f/' > f/. For U' < U the Kondo scale 
increases from the SU{2) value obtaining at U' = 0, to a much larger value at the 
SU{4:) point U' = U. For C/' > U — i.e. in the charge-Kondo regime — the scale then 
decreases very rapidly but continuously with increasing U' , vanishing at the critical 
point U' = U'^. We now consider each regime in turn. 

For U' = 0, the fJ-dependence of the SU{2) Kondo scale is of course well known, 
and our numerical results recover correctly the asymptotic strong coupling form [2, 29] 

Tr^^)ocrf/iexp(^^j. (6.1) 

On the SU{4) line U' = U hy contrast we find 

Tr^'ocrf/texpj^^j (6.2) 

for large U (including both the exponential and the prefactor), where the reduction 
in the exponential argument by a factor of 2 = ^ is as expected for SU {2N) Kondo 
behaviour [2] . The marked enhancement of Tk seen in figure 5 as f/' — > [/ thus reflects 
rpSU{4) ^ jj,S[/(2)ji/2 ij ^ I jjj fa^Q^ .^g ggg from the figure that Tk grows most 

rapidly in the close vicinity of the 5'f/(4) point; and for U — U' ':$> T^^^'' is not 
substantially difi'erent from its value at U' = 0. In later sections we shall see that this 
is also reflected in essential persistence of the U' = Q SU (2) spin-Kondo physics over 
a wide range of U' in the SC phase. 
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The second SC regime of interest indicated by figure 5 is the charge-Kondo regime 

U < U' < U'^, where the Kondo scale drops rapidly with increasing U' and ultimately 
vanishes at the phase boundary. The vanishing of Tr- is found to be described very 
accurately by the Kosterlitz-Thouless (KT) [36] form 

Tk = c{U)r exp I ^^EL I (6.3) 

where a and c are positive constants for a given U. We note too that (6.3) holds for 
all U > (and not just strong-coupling), whence the KT behaviour near the phase 
boundary is generic. The inset to figure 5 shows a typical fit to the KT form for U ~ 7, 
where the constants a and c, and the critical U^, have been chosen to achieve the best 
possible fit. It is evident from the figure that the numerical results agree very well with 
equation (6.3) over many orders of magnitude, and that the KT behaviour is in fact 
observed over much of the charge-Kondo regime. The specific value of [/^ ~ 7.046 thus 
found also agrees very well with that obtained from the bisection method discussed 
at the beginning of this section, indicating as such the rcliablity of the NRG results. 

In the CO phase by contrast there is no low-energy Kondo scale at all: as discussed 
in the following sections, the low-energy physics here is controlled simply hy U' — U , 
which remains finite when the critical is approached from above. 

6.2. Thermodynamics 

We now turn to a more detailed analysis of the DQD thermodynamics. This enables 
us both to verify the schematic NRG flow diagram shown in figure 2, and to determine 
the temperature scales at which the flows pass from one fixed point to the next. In the 
SC phase, it is of course Tk that sets the scale for the ultimate crossover to the SC 
fixed point. In strong coupling, the low-temperature physics in this phase can thus be 
described in terms of universal functions of T/Tk; we therefore end the section with 
analysis of this scaling behaviour. 

6.2.1. Entropy. Figure 6 shows S = S'lmp (see section 3.3) versus T/T for a range 
of different interdot interactions U' at fixed U = 7; we consider points in both the 
SC (U' < C/g ~ 7.046) and CO (U' > U'^) phases. In each case the entropy decreases 
stepwise with temperature, beginning at the high-T asymptote of In 16 corresponding 
to the FO fixed point [29]. The plateaus seen as the temperature is lowered (i.e. 
iteration number N increased, T ~ DA'^/'^ [28, 29]) arise when the NRG flow is 
close to the fixed points of the system, and by identifying the characteristic entropy 
associated with each fixed point it is possible to deduce the route taken in each case by 
the NRG flow. The reader may find it helpful to refer to the schematic flow diagram 
of figure 2 in combination with the following discussion. 

We first recap briefly how the entropy varies with temperature in the well- 
understood C/' = limit [29], shown in flgure 6(a). With decreasing temperature it 
changes from In 16 to In 4 (= 2 In 2) on the scale T ^ U, reflecting the crossover from 
the FO fixed point to the quadruply-degenerate LM'^^(^) where charge fluctuations 
on the dots are frozen-out. The entropy then decreases logarithmically slowly with T, 

SU(2) 

until at T ~ , there is a crossover to the SC fixed point where S'lmp = 0. 

How is this changed at finite U'7 Consider first the S't/(4) point U' = U, 
figure 6(d). Here there is again just a two-stage crossover in the entropy on reducing 
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Figure 6. Impurity entropy S{T) versus T for U = 7 with (a) U' = 0, (b) 
U' = 6 (SC phase, spiii-Kondo regime), (c) U' = 6.9, (d) U' = 7 (5(7(4)), (e) 
U' = 7.03 (SC, charge-Kondo regime) and (f) U' = 8 (CO phase). The critical 
~ 7.046. Dotted lines show lnl6, ln6, ln4 and ln2, associated with the FO, 
LM'^^W, LM^^^(2)a^jj(j QQ gxed points respectively. 



T, but this time the plateau at intermediate-T is at In 6 rather than the In 4 seen for 
U' = 0; on these temperature scales the physics is controlled by the LM'^^^'*) fixed 
point (see figure 2), where all four states of the (nL,nij) = (1,1) dot configuration, 
and the pair of states (2,0) and (0,2), arc degenerate. The subsequent crossover from 
LM^^Wto SC occurs on a much higher temperature scale than the analogous crossover 
from LM'^'^(^) to SC in figure 6(a), consistent with the marked increase in the SU{4:) 
Kondo scale t^'^^^ over that oiU' — Q. 

If one now moves slightly away from U' = U, the situation becomes more 
complicated. On temperature scales T » \U' — U\, the form of the entropy is found to 
be the same as at the SU (4) point itself, but now on lower temperature scales there 
are additional crossovers to the other unstable fixed points of the system. For U' < U, 
e.g. U' = 6.9 as shown in figure 6(c), the NRG flow crosses from LM^^W lm^^(2) 
on the scale T ^ U — U', and hence the entropy drops from In 6 to In 4. (At lower 
temperatures still, the local moments of LM'^^^^^ are Kondo quenched and the entropy 
crosses over to the SC value of zero on the scale T ~ Tk as before.) On the other 
hand, if U' > U, it is the CO fixed point that is approached after LM'^'^^'*), on the 
scale T ^ U' — U. Here, the remaining dot degree-of-freedom is the charge psuedospin 
made up of the (2,0) and (0,2) configurations, and the entropy thus plateaus at In 2. 
What then happens at lower temperatures depends on the size of U' . For U' < t/^', as 
seen in figure 6(e) for U' = 7.03 (C/^ ~ 7.046), the CO fixed point is unstable and the 
flow eventually goes to SC (reflecting the quenching of the charge psuedospin on the 
temperature scale T ~ Tk (6.3)). For U' > by contrast, as shown in flgure 6(f), 



Renormalization group study of capacitively coupled double quantum dots 



20 



the CO fixed point is stable and its In 2 entropy thus persists down to T = 0. 

6.2.2. Susceptibilities We consider now the impurity spin- and staggered charge- 
susceptibilities defined at the end of section 3.3. These probe different aspects of the 
dot degrees-of-freedom, and show clearly that the transition is from the charge-Kondo 
state below to the doubly-degenerate charge ordered state above U^, the dot spin 
degrees-of-freedom playing essentially no role. 

Figure 7 shows results for the spin susceptibility. The main figure shows Txs 
versus T for a fixed U ~ 7 and a range of U' spanning the transition (U^ ~ 7.046): 
U' = 6, 6.9, 6.97, 7, 7.03, 7.1 and 8 (from top to bottom). In all cases the behaviour 
at high temperatures (T ^ U) is clear: each dot contributes 'half a Curie-law' to the 
spin-susceptibility, hence Txs ~ j- For f/' = 6, on lowering the temperature the form 
of the susceptibility is seen to be essentially a doubled version of that found for the 
the U' = limit [29], the U' = physics thus persisting far into the SC phase as seen 
in the previous section (c/ e.g. figure 6(a) and (b)). The susceptibility tends to the 
full Curie-law value of | as the NRG flow approaches LM'^'^^^^ , and then drops to zero 
at T ~ Tk where the free spins are quenched by the Kondo effect. 

As U' is increased toward U, Xs begins to show signs of the LM'^^^^^ fixed point; 
this is readily shown to have a characteristic Txs = |) and is seen as a shoulder in 
the U' = 6.9 curve as it heads toward to the LM'5''^(2) i^^^^ i_ U' = U itself, 
this shoulder disappears and Txs evolves smoothly to zero with decreasing T. 

Perhaps the most striking feature of figure 7 is that it shows no sign of the 
quantum phase transition: the spin susceptibilities evolve smoothly as U' passes 
through U'-. and as T ^ 0, Txs vanishes in both phases. The reason is as follows. 
For values of U' much more than ~ rpSUii) ^^^ow U, the low-temperature physics 
consists of the quenching of the LM'^'^f^^ local moments on the temperature scale 
Tk ~ T^''^^\ On increasing U' towards the SU{A) point U' — U, all that happens is 
that Tk rises slowly (see figure 5) and hence Txs crosses over to zero at progressively 
higher temperatures - as seen in figure 7. Once U' reaches U, it is the quenching of 
the 'magnetic' {nL,nii) = (1, 1) components of the LM''''^'^^^ local-moment states that 
show up in the spm-susceptibility, as a crossover from i to at T ^ j,su{i) ^ jj, 
now increased further, into the charge-Kondo regime, the magnetic degrees-of-freedom 
are frozen out at an even higher scale T U' — U (the difference in energy between 
the (1, 1) and (2,0)/(0, 2) dot states), corresponding to the crossover from LM'^^^^'Ho 
CO. Whether the flow eventually heads to SC at very low temperatures or remains 
at CO is irrelevant: the magnetic behaviour is already lost at T ^ f/' — J7, and hence 
the susceptibility evolves smoothly through U' = U^. 

This is further confirmed by the inset to figure 7, showing the zero-temperature 
XsiT = 0) versus U'. We find that Xs(0) crosses over from a relatively large value 
indicative of the SU (2) x SU (2) spin-Kondo physics to a somewhat smaller value at 
the SU (4) point, and then decreases further still as the charge-Kondo regime is entered 
and the non-magnetic (2, 0) and (0, 2) impurity configurations dominate the T ^ 
behaviour. Upon moving into the CO phase, the spin-susceptibility simply evolves 
continuously, and is found to decay smoothly to zero as ?7' — > oo. 

The situation is quite different for the staggered charge susceptibility Xc ((311))) 
the charge pseudospin analogue of the spin susceptibility. Figure 8 (left) shows 
Tx7{T) versus T/T for fixed U = 7, with U' = 7, 7.04, 7.044, 7.048, 7.1 and 
8. Here, the behaviour characteristic of the two phases is qualitatively distinct. 
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Figure 7. Spin susceptibility Txs{T) versus T/r for [/ = 7 and U' = 6, 6.9, 
6.97, 7, 7.03, 7.1 and 8 (from top to bottom). Inset: the zero-temperature Dxs(O) 
versus U', again for [7 = 7. 
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Figure 8. Left: staggered charge susceptibihty. Txc (T) versus T/F, for fixed 
(7 = 7, with U' = 7, 7.04, 7.044, 7.048, 7.1 and 8 (from bottom to top) and 
~ 7.046. Right: [/'-dependence of the Wilson ratios iij (sohd Une) and Rs 
(dashed) for U = 7. 

with TXc vanishing as T ^ in the SC phase but tending to a finite constant 
in the CO phase. This reflects the fact that the charge pseudospin comprised of 
the (2,0) and (0,2) states is screened by the charge-Kondo eSect at low-T in the 
SC phase, but remains free in the CO phase and hence gives rise to a staggered 
charge susceptibihty oc So whereas the spin-susceptibihty Xs(0) shown in 

the inset of figure 7 evolves smoothly through the QPT, the T = staggered charge- 
susceptibility diverges as U' U'~ , and remains infinite throughout the CO phase. 
This is directly analogous to the magnetic susceptibility of the two-dimensional xy 
model [36], the prototype of Kosterlitz-Thouless physics. Note moreover that the 
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divergence of Xc (0) tracks precisely the vanishing of the Kondo scale as the transition 
is approached: from (5.9) with = I/7, and using R~ — > 4 as f' ^ ((5.15)), 
we have Xci^ = 0) ~ S/^tt^T^). 

The full T-dependences of the staggered charge susceptibilities shown in figure 8 
can likewise be explained in terms of the fixed-point picture used for Xs- Here we 
simply note two points. First, the crossover in TxJ to zero occurs on the scale T ~ Tk 
only when U' > U; for U' < U the charge degrees of freedom are frozen-out on the 
scale T ^ U — U' instead (cf. the discussion of the spin-susceptibility above). Second, 
the limiting values of TXc as T ^ in the CO phase are clearly [/'-dependent; this 
arises naturally from the potential scattering generated at the CO fixed point (and 
considered further in the Appendix) . In the absence of such potential scattering — as 
arises in the limit £/' ^ 00 — the T ^ limit of Tx~ is unity. 

6.2.3. Wilson Ratios The three Wilson ratios characteristic of the SC phase {Rg and 
R^) have been considered in section 5.2, and exact results for them obtained at the 
points U' = 0, U and U^~. Here we show the corresponding numerical results as a 
function of U'; this shows the U' scales over which the ratios change from one value 
to the next, and also provides an indication of the accuracy of the NRG calculations. 
Figure 8 (right) shows the [/'-dependence of Rg and R~ for U = 7; as U here is large, 
the uniform charge susceptibility and hence is found in practice to be essentially 
zero over the entire U' range. We see that the Rg = 2, R~ = result known for the 
[/' = limit [28, 37] persists with increasing U' until very close to the SU{A) point; 
but then the two Wilson ratios undergo a rapid crossover to the values i?,, = and 
R~ = 4 ((5.15)), which values obtain in practice over a large part of the charge-Kondo 
regime. At the SU (4) point itself, we recover numerically (to within ~ 2%) the exact 
result Rs = R~ = | ((5.13)) expected from the symmetry of the model. 

6.2.4. Universal scaling behaviour. Finally, we consider the important issue of the 
scaling properties of thermodynamics in the SC phase (naturally in strong coupling, 
/7 3> 1, where the Kondo scale Tk is exponentially small and well separated from 
all non- universal scales). We point out immediately that it is not possible to define 
a single universal scaling function (for each thermodynamic property) that applies 
throughout the SC phase. The reason is simple: although the spin-Kondo regime 
([/' < U), SU{A) point ([/' = U) and charge-Kondo regime ([/<[/'< U'j all belong 
to the SC phase, the crossover from one to the next being continuous with increasing 
U' , each has its own characteristic low-temperature behaviour because the particular 
degrees-of-freedom quenched by the Kondo effect are different in each case. This is 
particularly apparent from the Wilson ratios shown in figure 8. We thus focus on the 
three distinct regimes separately, beginning with spin-Kondo. 

In the spin-Kondo regime, it is the crossover of the spin-susceptibility Txs on the 
scale of Tk that characterises the physics of the spin-Kondo effect. In figure 9, we 
show Txs versus T/Tk for U' = Q (solid line), 6 (dotted line) and 6.9 (dashed line), 
all with U = 7. It is clear from the figure that the rescaled susceptibilities all collapse 
onto a common curve at low temperatures (and that the crossover from LM'^^*^^) to SC 
indeed occurs at T ~ Tk)- Also plotted in figure 9 (as a set of points) are double the 
universal sm^^e-impurity Kondo susceptibilities taken from Table V of [29] , which are 
clearly in excellent agreement with the scaling curves obtained here. [The Kondo scale 
of [29], here denoted by T\y, is easily shown to be related to our definition Tk = 1/7 
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Figure 9. Universal scaling of the spin-susceptibility in the spin-Kondo regime: 
Txe versus T/Tk for [/ = 7 and U' =0 (solid hne), 6 (dotted line) and 6.9 (dashed 
line). The points plotted in the figure are double the universal single-impurity 
Kondo susceptibilities, taljen from Table V of [29]. 

by Tk = 3Tw/{tt'^w), with w = 0.41071 ... the Wilson number [2]]. 

In the charge-Kondo regime, it is by contrast the staggered charge-susceptibility 
that shows universality. In figure 10 we plot Tx~ versus T/Tk for U = 7 with 
U' = 7.044 (sohd hne), 7.04 (dotted) and 7.03 (dashed), and also for [/ = 8 with 
U' = 8.028 (dot-dashed line). That the curves collapse onto a common form at low- 
temperatures is again evident, the scaling becoming better as U' U'r where the 
fixed point Hamiltonian (5.1) becomes increasingly valid. Points in the spin-Kondo 
regime {W < U) do not scale onto this universal curve, because there the charge 
susceptibility vanishes on a finite energy scale U — U' (such that {U — U')/Tk — > oo 
is thus 'projected out' in the formal scaling limit Tk 0). 

Finally, the behaviour along the S'J7(4) line U' = U is seen in figure 11 where 
Txs versus T/Tk is shown for U = tj' = 7 (solid hne), 8 (dotted) and 9 (dashed); 
and from which universality is again evident. By exploiting the SU (4) symmetry the 
behaviour of versus T/Tk is identical here, and that is indeed recovered by the 
numerical results. 

7. Concluding remarks 

Using the numerical renormalization group we have considered a symmetrical, capac- 
itively coupled DQD in the two-electron regime, and investigated the evolution of the 
system as a function of both the interdot and intradot coupling strengths. The range 
of physical behaviour found is both broad and subtle, encompassing a strong coupling, 
Fermi liquid phase - itself divided into spin-Kondo and charge-Kondo regimes sepa- 
rated by a spin-charge entangled SU{4:) line - as well as a broken symmetry charge 
ordered phase, and the Kosterlitz-Thouless quantum phase transition occurring be- 
tween the two. As far as specific physical properties are concerned, we have focussed 
largely in this paper on static properties, such as 'impurity' thermodynamics, spin- 



Renormalization group study of capacitively coupled double quantum dots 



24 




10 10 10 10 



10 10 10 10 10 10 



T/Tk 



Figure 10. Universal scaling of the staggered charge-susceptibility in the charge- 
Kondo regime: Tx7 versus T/Tk for [/ = 7 with U' = 7.044 (soUd line), 7.04 
(dotted) and 7.03 (dashed), and for {/ = 8 with U' = 8.028 (dot-dashed Une). 




10 10 10 10 10 10 



T/Tk 

Figure 11. Universal scaling of the spin-susccptibility along the SU(A) line: Txs 
versus T/Tk for U = U' = 7 (solid line), 8 (dotted line) and 9 (dashed line). In 
this case the spin and staggered charge susceptibilities coincide, Xs = Xc ■ 



and charge-susceptibilities etc. Importantly, the NRG approach also enables us to 

determine dv d properties such as local single-particle spectra [39] and related 

transport properties. It is to these we shall turn in a subsequent paper [31]. 
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Appendix 

In the strongly correlated regime, one can of course obtain a low-energy effective model 
for the n = 2 regime by performing a Schrieffer- Wolff transformation [2] of (2. la, 6) to 
order V'^. The final result is 



where 



i.k.k ,(T 



+ aJ + K) V 4,q,aLa,^,,, (A.2) 



i,j^i,k,k' ,a,a' 



and the coupling constants are 



4y2 oy^ 
AcP = U'-U, J=— and if = __. (A.3) 

The first term on the right-hand side of (A.2) acts solely on states of the (ni,n/{) = 
(1,1) dot configuration. It is simply a separable pair of spin-i Kondo couplings 
between each dot and its respective lead. By contrast, the second term in (A.2) acts 
only on states of the (0, 2) and (2, 0) configurations; this is the correlated potential 
scattering term that arises at the CO fixed point (see section 4.1). It is the third 
term in (A.2) that mixes the different charge sectors, but note that although (1,1) 
connects to both (0, 2) and (2, 0) to order V^, the (0, 2) and (2, 0) states do not connect 
directly to each other. For the special case of U' = U, one obtains K = and the 
Hamiltonian then takes the manifestly SU (4) symmetric form 

k,m k,k' ,m,m' 

with m = as before. 

The effective Hamiltonian above leads to a simple physical picture of the DQD 
physics with increasing U' . When U' = 0, only the first term in (A.2) need be 
retained: the (0,2)/(2,0) states lie an energy U above the (1,1) ground state by 
virtue of the l\(p term in (A.l), and hence do not contribute to leading order in 
y^/?7. The effective Hamiltonian is thus two uncoupled spin-i Kondo models as 
one would of course expect. This situation does not change with increasing U' until 
the energy gained by mixing (0, 2)/(2,0) into the (1, 1) ground state — of order the 
SU (4) Kondo temperature j"^^^^ - jg sufficient to outweigh the cost of promoting 
electrons into these excited states, A(^. Hence one would expect a rapid crossover 
from the SU{2) x SU{2) Kondo effect characteristic of U' = 0, to the 5C/(4) Kondo 
behaviour of [/' = [/, at a J7' ^ t/ — 0{T^^'"^^) as indeed seen in figure 5. Once 
U' > U, the order of the (0, 2)/(2, 0) and (1,1) states is reversed with the latter 
being the excited state through which the former interconvert (via the third term in 
(A.2)). If A(f) is sufficiently small, then the intcrconvcrsion of (0, 2) and (2, 0) via (1, 1) 
gives rise to a charge-Kondo effect that quenches the (2, 0)/(0,2) charge pseudospin, 
and thus the SC phase persists above U' (as shown in figure 4 and the associated 
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discussion). Once U' > U + 0(T^^^^^) however, the stabilisation energy associated 
with this chargc-Kondo cfFcct cannot compensate for the additional cost of occupying 
the (1, 1) excited states, the (0, 2) and (2, 0) states do not therefore communicate, and 
the system undergoes the quantum phase transition to the CO phase. 

Deep in the CO phase, when U' 3> U, one can simphfy the effective Hamiltonian 
further still by neglecting the terms in (A. 2) that involve the (1, 1) states. This leaves 
only the correlated potential scattering described by the second term. The resulting 
effective Hamiltonian in the (0, 2)/(2,0) manifold is then a contimmm form of the 
CO fixed point Hamiltonian (4.3), as mentioned in section 4.1. This means that one 
can obtain the fixed point Hamiltonian (4.3) by direct discretisation of the effective 
Hamiltonian; from which the following relationship between K and K is obtained 

K = AapK (A.5) 

with Af^ = 5(1 + A~^)(1 — A~^)~^ln(A) the standard correction factor for the 
discretisation of the conduction band [29]. In terms of the parameters of the full 
DQD Hamiltonian (2.1a,Z)), this becomes (using (A. 3)) 

r> 41nA 1 

K = =r, (A.6 

7r2(l-A-i) 2U'-U 

which we plot in figure Al (dashed line) as a function of U' for fixed U = 7 in the 
CO phase. The solid line shows the value of K obtained from our full numerics, and 
is indeed seen to fall rapidly onto (A.6) as one moves deeper into the CO phase. The 
inset of figure Al shows the corresponding anomalous exponent —a, calculated from 
K using (4.4) and (4.6). At the phase boimdary itself (t/^ ~ 7.046), a vanishes and the 
perturbation dHio in table 2 is marginal. On moving into the CO phase, a increases, 
tending to | in the limit U' ^ 00. Note that K can of course be obtained by analysing 
the NRG flows directly; in practice the result is indistinguishable from that obtained 
using the method described above. 
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Figure Al. CO fixed point potential scattering strength, K, vs. U' for fixed 
U = 7. The dotted line marks the critical C/^ ~ 7.046 while the dashed line is the 
approximate result (A.6) derived from the 0{V'^) eflfective Hamiltonian. Inset: 
the corresponding anomalous exponent —a, calculated using (4.4) and (4.6). 
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